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We present empirical relations for determining the amount by which the effective temperatures and radii — and therefore 
the estimated masses — of low-mass stars and brown dwarfs are altered due to chromospheric activity. Our relations are 
based on a large set of well studied low-mass stars in the field and on a set of benchmark low-mass eclipsing binaries. 
The relations link the amount by which an active object's temperature is suppressed, and its radius inflated, to the strength 
of its Ha emission. These relations are found to approximately preserve bolometric luminosity. We apply these relations 
to the peculiar brown-dwarf eclipsing binary 2M0535— 05, in which the active, higher-mass brown dwarf has a cooler 
temperature than its inactive, lower-mass companion. The relations correctly reproduce the observed temperatures and 
radii of 2M0535— 05 after accounting for the Ha emission; 2M0535— 05 would be in precise agreement with theoretical 
isochrones were it inactive. The relations that we present are applicable to brown dwarfs and low-mass stars with masses 
below 0.8 Mq and for which the activity, as measured by Ha, is in the range —4.6 < log Ln a /Lboi < —3.3. We expect 
these relations to be most useful for correcting radius and mass estimates of low-mass stars and brown dwarfs over their 
active lifetimes (few Gyr). We also discuss the implications of this work for determinations of young cluster IMFs. 
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1 Introduction 

Observational evidence strongly indicates that the funda- 
mental properties of low-mass stars can be altered in the 
presence of strong magnetic activity (e.g., Morales et al. 
2008; Lopez-Morales 2007; Ribas 2006). In particular, ob- 
servations of active, low-mass eclipsing binary (EB) stars 
have found the empirically measured stellar radii (R) to be 
inflated by «10%, and the empirically measured stellar ef- 
fective temperatures (T ff) to be suppressed by «5%, rel- 
ative to the predictions of standard theoretical stellar evo- 
lution models, which better match the properties of inac- 
tive objects (see Coughlin et al. 2011; Kraus et al. 2011; 
Morales et al. 2010; Stassun et al. 2009, and references 
therein). Resolving these discrepancies will be critical to the 
ongoing development of accurate theoretical stellar models 
(see Stassun et al. 2010). Accurate estimates of stellar radii 
are especially important in the context of searches for tran- 
siting exoplanets, which rely upon the assumed stellar ra- 
dius/density to infer the planet radius/density. 

Activity effects also lead to errors in object masses (M) 
when these are derived from T e g . A particularly salient ex- 
ample is 2M0535-05, an EB in the Orion Nebula Cluster 
(age ~1 Myr) comprising two brown dwarfs (Stassun et 
al. 2006, 2007). The primary and secondary brown dwarf 
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(BD) components of 2M0535— 05 have dynamically mea- 
sured masses of 60±3 and 39±2 Mj up , respectively, and 
T off ratio of T ± /T 2 = 0.952 ±0.004 (Gomez Maqueo Chew 
et al. 2009). That is, the system exhibits a reversal of the 
usual M-T e ff relation, such that the primary component is 
cooler than its companion. Figure[T]shows the 2M0535— 05 
system on the Hertzsprung-Russell (H-R) diagram. The sec- 
ondary BD's T e ff and bolometric luminosity (L\, a i, calcu- 
lated directly from the empirically measured T c & and R) 
place it at a position that is consistent with that predicted 
by the model isochrone. In contrast, the primary is far dis- 
placed from its expected position, and so appears to have 
a mass of only ~25 Mj up — more than a factor of 2 lower 
than its true mass — on the basis of its low T c & . Reiners et al. 
(2007) used spectrally resolved Ha measurements to show 
that, whereas the secondary BD in 2M0535— 05 is chromo- 
spherically quiet, the primary BD is highly chromospheri- 
cally active, perhaps a consequence of its rapid rotation. 

Since magnetic activity seems to alter the fundamen- 
tal properties of both stars and BDs, it would be valuable 
to have an easily observable metric with which to quan- 
titatively assess the degree to which a given object's T e ff 
has been suppressed and its radius inflated. Here we present 
such an empirical metric by relating the degree of T c ff sup- 
pression and radius inflation to the strength of Ha emission, 
a commonly used and readily observable tracer of chromo- 
spheric activity (Scholz et al. 2007; Berger 2006). 
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Fig. 1 H-R diagram for the primary and secondary BD 
components of the EB 2M0535-05 in the —l-Myr Orion 
Nebula Cluster (Stassun et al. 2006, 2007). The measured 
T e g and Lboi for both BDs are represented as blue sym- 
bols. Measurement uncertainties in T c s and Lboi are rep- 
resented by the error bars at upper right. The dynamically 
measured masses of the primary and secondary are repre- 
sented as red bars on the 1-Myr theoretical isochrone of 
Baraffe et al. (1998). The measured Lna/iboi for the two 
components are indicated next to the blue symbols. The ac- 
tive primary appears far cooler than expected, and therefore 
appears to be much younger than the secondary and to have 
a mass of only w27 Mj up based on its observed T c g , a fac- 
tor of 2 lower than its true mass. Shifting the position of the 
active primary (arrow) using our empirically calibrated Ha- 
based relations brings the primary into much closer agree- 
ment with its theoretically expected position in the HR dia- 
gram (black symbol); this is where the active primary would 
be if it were not active. 



2 Methods and Data Used 

We use two different samples and approaches to empirically 
determine a relationship between T c g suppression, R infla- 
tion, and the level of activity as measured from Ha. 

First we use the large set of nearby field M dwarfs 
with well measured spectral types and Ha equivalent widths 
(EWs) from the PMSU catalog (Reid et al. 1995; Hawley et 
al. 1996). Following Morales et al. (2008), we restrict our- 
selves to the sample of 746 stars with distances determined 
directly from trigonometric parallaxes. Fig. |2] shows the es- 
timated T c ff and R of the PMSU sample stars as a function 
of M. Strongly chromospherically active stars — defined as 
those with Ha in emission — show a clear displacement to 
lower T ff and to larger R relative to both the theoretical 
isochrone and to the non-active stars, whereas non-active 
stars more closely track the isochrone. From these T c ff and 



R offsets, and the observed Ha emission, we derive linear 
relationships between AT c ff, AR, and log Lna/^boi (see 
Fig. [3]). The relationships are statistically significant with 
>90% confidence. 
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Fig. 2 T ff vs. Mass (top) and Radius vs. Mass (bot- 
tom) for M-dwarfs with trigonometric distances and Ha 
measurements from the PMSU catalog. Active objects (Ha 
in emission) are filled symbols. Solid curve is a 3 Gyr 
isochrone (Baraffe et al. 1998). Dashed curve is a polyno- 
mial fit to the non-active objects. Ha-active dwarfs are sig- 
nificantly displaced to lower T c s and larger radii compared 
to both the isochrone and to the non-active dwarfs. 



The second approach uses the much smaller sample 
of stars in EBs that have directly measured masses, radii, 
and reliable T ft, but for which we must use X-ray flux 
as a proxy for Ha. We use the small set of low-mass EBs 
with accurately measured Al, R, T c g, and X-ray luminosi- 
ties (Lx) from Lopez-Morales (2007). The sample includes 
11 individual stars in 7 EB systems spanning the range 
T cS =3 125-5300 K and A/=0.21-0.96 M . We begin with 
the correlation of AR vs. LxlLho\ already demonstrated in 
that work, which we rederived using the fundamental stellar 
data compiled in Lopez-Morales (2007) and the same 3 Gyr 
isochrone of Baraffe et al. (1998) as above. Lopez-Morales 
(2007) did not discuss the complementary correlation with 
AT c ff, but this information is also contained in the EB data, 
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so we also derive the relationship AT e g vs. LxlLhoh Next 
we use the empirical Lx/L^ \ vs. L^a/L^d relationships of 
Scholz et al. (2007) and Delfosse et al. (1998) convert the 
above from relationships on L^/Lboi into relationships on 
Ln a /L\> \. The resulting relationships between AT c ff , AR, 
and log ina/iboi (Fig. [3j> are statistically significant with 
>95% confidence. 




'°9 L Ha /L bol 

Fig. 3 T e g suppression (top) and radius inflation (bot- 
tom) as a function of fractional Ha luminosity for both field 
M-dwarfs (filled symbols) and low-mass eclipsing binaries 
(error bars). Linear fit relations to both sets are dashed lines. 
The final averaged best-fit relation the solid line. See text for 
the linear fit coefficients. 



3 Empirical Relations Linking Radius 
Inflation and T cff Suppression to Ha Emission 

Combining the results from the field M dwarfs and EBs, we 
obtain the following final relations (in percent units): 

AX 1 

^ = (-4.71±2.33)x(log W£boi+4)+(-4.4±0.6)l 

s 



AR 

As an example of how the apparent mass of an object 
can be altered by activity, we calculate how the observed 



(15.37±2.91)x (logL HQ /iboi + 4) + (7.1±0.6). 



H-R diagram positions of the two BDs in 2M0535— 05 are 
altered by the above relations. In effect, we are seeing how 
the 2M0535— 05 system would appear in the H-R diagram 
were the system completely inactive. We use the Ln a /L^ \ 
measurements of Reiners et al. (2007) shown in Fig.Q] We 
find that the primary BD has been displaced by AT e ff = 
-6.9 ± 1.4% and AR = 15.2 ± 1.7%, and the secondary 
BD by at most AT off = -3.0±0.9%and AR < 2.5±1.1% 
(implying nearly constant Lboi for both). Shifting the posi- 
tion of the primary BD accordingly (Fig.Q]) shows that, if it 
were inactive, it would be in excellent agreement with the 
theoretically expected position for its known mass. 

4 Impact on Inferred Initial Mass Functions 

The vast majority of masses for stars and BDs in young clus- 
ters can only be determined by comparison with theoretical 
evolutionary tracks in H-R diagrams, using either Lboi or 
T e g or both (e.g., Scholz et al. 2012). Therefore, if there 
is indeed a relation between magnetic activity and R in- 
flation / T c ff suppression, this will also affect estimates of 
stellar and substellar masses derived from T G ff , especially at 
young ages when activity levels are high. Here we examine 
the two ways of estimating masses — T ff and Lboi — that are 
commonly used in the literature and investigate the impact 
of magnetic activity on the derived masses. 

First we use T c s to estimate masses from model 
isochrones at 1 Myr. At each model mass, we consider a 
range of activity levels and apply offsets to the model T c g 
based on our empirical relations above. We then use these 
suppressed T c ft values to estimate the masses that would 
be inferred from the isochrone. As expected, this proce- 
dure leads to a systematic underestimation of the masses. 
At high levels of activity, the effect can be substantial. For 
log L-aJL^oi = —3.3, which corresponds to the saturation 
limit in young low-mass stars and in young associations 
(Scholz et al. 2007), the mass estimates are a factor of ~ 2 
lower than for objects with low levels of magnetic activity 
(log LHa/Lbol < —4.5). Note that because the mass-T e ff 
relationship is less steep at older ages, the underestimation 
of mass is less severe at older ages (> 100 Myr). 

In the second test case, we use Lboi derived from K- 
band absolute magnitudes to derive masses. In this case the 
influence of magnetic activity is introduced by the T c g de- 
pendence of the bolometric correction, since Lboi is (as dis- 
cussed above) practically not affected. We compute spectral 
types from the suppressed values of T e g, and then compute 
K -band bolometric corrections from the spectral types, and 
we combine these with the model A'-band absolute magni- 
tudes to find Lboi- Finally, we estimated masses from the 
model isochrone and Lboi- This second method still leads 
to an underestimate of the masses, but the effect is much 
smaller than when directly estimating masses from T e s. The 
change of the bolometric corrections with increasing mag- 
netic activity is quite small (at most 10%), resulting in rela- 
tively minor changes of a few percent in the mass estimate. 
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An important application of T c g -based mass estimates is 
the determination of initial mass functions (IMFs) for young 
clusters. M-dwarfs with ages from 10 to 100 Myr exhibit Ha 
emission in the range log Ln a /Lboi = —4.2 to —3.3 (Scholz 
et al. 2007). Thus, our analysis implies that objects inferred 
to be BDs from T c ff-based (or spectral type based) mass 
estimates could actually be low-mass stars since masses will 
be underestimated by up to a factor of 2. 

We further evaluated the impact of this effect on mea- 
surements of the slope of the IMF, a (in dN/ dM oc M ~ Q ). 
For this purpose, we assume a measured slope of a = 0.6 
(see Scholz et al. 2012), calculate the IMF based on that 
slope for M < 0.6 M©, correct the masses for a given 
level of Ha emission using our relations above, and re- 
determine the slope a. If we start with log L^ a IL^ \ = —4.0 
for BDs and —3.5 for low-mass stars (Barrado et al. 2003), 
the masses have to be corrected by factors of 1.3-1.7 for 
BDs, and by 2.2-2.5 for low-mass stars. As a result, the 
slope changes from a = 0.6 to a = 0.5. In general, we ex- 
pect that a will be underestimated by <; 0.1, if the masses 
are estimated from T e s and activity is not taken into ac- 
count. In addition, the peak mass of the IMF would be un- 
derestimated by up to a factor of 2. 

5 Discussion 

It remains unclear what is the underlying physical mecha- 
nism driving the correlations between AT e s, Ai?, and Lho- 
Chabrier et al. (2007) and MacDonald & Mullan (2009) 
have suggested that a sufficiently strong field could suppress 
convection, inhibit heat transfer, and thus inflate (and cool) 
the stellar surface. Since such a field would also likely re- 
sult in chromospheric activity, one might therefore expect 
the correlations with Ha emission that we have derived. 

Browning (2008) has shown in global numerical models 
that fully convective stars can host large kG strength fields. 
However, such fields alone appear to be too weak to pro- 
duce the radius inflation and T e s suppression observed in 
2M0535-05 (MacDonald & Mullan 2009). Chabrier et al. 
(2007) also suggest that fully convective objects should be 
less affected by the same convective inefficiencies invoked 
to explain radius discrepancies at higher masses. It is pos- 
sible that a combination of rotation and magnetic activity 
contribute to both inflation/suppression and Ha emission in 
such a way as to produce our empirical relation without a 
causal correspondence between Ha and the magnetic field. 

An alternative explanation for T c s suppression is a spot 
covered surface. Since spot coverage is also controlled by 
magnetic fields, a correlation with Ha might still exist. 
However, Mohanty et al. (2010) and Mohanty & Stassun 
(2012) have now shown from a spectral fine-analysis of the 
2M0535— 05 system, observed at high resolving power dur- 
ing both primary and secondary eclipses, that such a large- 
spot scenario is strongly disfavored as an explanation for the 
T e g suppression of the primary BD in 2M0535— 05. Thus, 



while the T c h suppression mechanism produces a clear cor- 
relation with chromospheric Ha activity as we have shown 
here, it evidently does not in all cases effect this correlation 
directly through surface spots. 

6 Conclusions 

We have shown that there exists a correlation between the 
strength of Ha emission in active M-dwarfs, and the de- 
gree to which their T e g are suppressed and radii inflated 
compared with inactive objects and theoretical evolutionary 
models. Our relationships above should prove directly use- 
ful for a variety of applications in which accurate estimates 
of masses and radii are needed for low-mass stars and BDs, 
and should assist in correcting inferred initial mass func- 
tions at young ages. 

While promising, the correlations we have derived con- 
tain significant scatter, and they are currently limited by the 
lack of a single sample of stars with both Ha and direct ra- 
dius measurements. We therefore encourage researchers to 
publish Ha measurements, as this is usually available from 
the spectra used to determine EB radial velocities. 

Finally, the relations we have determined already in- 
dicate quite clearly that the radius inflation and tempera- 
ture suppression mechanism operates in such a way that the 
temperature suppression and radius inflation almost exactly 
cancel in terms of their effect on the bolometric luminosity. 
Moreover, the relations between activity, T e g suppression, 
and radius inflation do not appear to manifest any obvious 
discontinuity across the fully convective transition (see also 
Stassun et al. 2010). These are important, fundamental clues 
to the physical nature of these effects, and should help to 
constrain theoretical models that are being developed to ex- 
plain these phenomena (e.g., Chabrier et al. 2007; MacDon- 
ald & Mullan 2009). 
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